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We study a model of frustration of decoherence in an open quantum system. Contrary to other 
dissipative ohmic impurity models, such as the Kondo model or the dissipative two-level system, the 
impurity model discussed here never presents overdamped dynamics even for strong coupling to the 
environment. We show that this unusual effect has its origins in the quantum mechanical nature 
of the coupling between the quantum impurity and the environment. We study the problem using 
analytic and numerical renormalization group methods and obtain expressions for the frequency and 
temperature dependence of the impurity susceptibility in different regimes. 
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I. INTRODUCTION 

In physics there is a large class of problems that can 
be described in terms of a single quantum mechanical de- 
gree of freedom interacting with an environment. Exam- 
ples range from magnetic impurities in metals, supercon- 
ductors, and magnets, macroscopic quantum tunneling 
in superconducting interference devices (SQUIDS) and 
molecular magnets^, to qubits in quantum computers^. 
The common thread between all these problems is the 
dramatic effect that the dissipation has on the quantum 
dynamics of the impurity 3 -. In particular, one of the most 
important effects of an environment on a quantum sys- 
tem is decoherence, that is, the destruction of quantum 
mechanical effects. Decoherence is the unavoidable con- 
sequence of the fact that no system in nature is really 
isolated. 

Impurity problems can be often reduced to an effec- 
tive one-dimensional boundary problem that allows the 
use of powerful non-perturbative theoretical techniques. 
The Kondo model is probably one of the best known im- 
purity problems and has been studied with a large num- 
ber of theoretical tools, from the exact solution via Bethe 
ansata^, numerical renormalization group£ , to conformal 
field theory^. The Kondo problem represents a universal- 
ity class of open quantum systems where dissipation and 
decoherence play a fundamental role. In its anisotropic 
form, the Kondo effect can be mapped via dimensional re- 
duction and abelian bosonization to the ohmic dissipative 
two-level system (DTLS) problem^. The Kondo effect 
can be thought as a situation where decoherence is ex- 
treme, in the sense that the spin is completely screened by 
the environmental excitations in the formation of the so- 
called Kondo singlet. Moreover, impurities can be used 
as probes for the understanding of the environment itself 
and in some cases can even determine the properties of 
the environment in a self-consistent manner. This occurs 
in the case of the dynamical mean-field theories (DMFT) 
where the solution of a many-body problem reduces to 



the solution of a self-consistent impurity problem^. Fur- 
thermore, systems where the competition between differ- 
ent phases of matter lead to the appearance of magnetic 
inhomogeneities (such as in the case of Griffiths-McCoy 
singularities in heavy fermion alloys) can many times be 
reduced to effective impurity problems^. 

In this paper we are going to describe a model for open 
quantum systems that cannot be described within the 
Kondo universality class. This model describes an effect 
that we call frustration of decoherence where decoher- 
ence is reduced by a pure quantum mechanical effect. It 
is important, therefore, that one understands the physics 
behind the standard model of dissipation described by 
the Kondo or the DTLS and how it relates to the prob- 
lem of decoherence. Since the connection between the 
Kondo problem and decoherence is not commonly dis- 
cussed in the literature we will review some of the key 
features of the DTLS and its connection with the prob- 
lem of decoherence. 

The DTLS can be described as a single spin half, S = 
(Si, S2, S3), coupled to a set of independent harmonic 
oscillators via the Hamiltonian (we use units such that 
h = l = fcs): 



-Hdtls = AS 3 + ~y=Si Vfc(afc - Qfc) 

V 2L fe> 

+ ^^vka\a k , (1) 

k 

where A is the tunnel splitting between the eigenvalues 
of Si, A is the coupling to an environment of bosons 
with one-dimensional momentum k, and energy disper- 
sion uik = vk (v is the velocity of the excitations that 
we set to unity, v = 1, from now on) and creation and 
annihilation operators a\ and a k , respectively (L is the 
linear size of the system). The operators obey canonical 
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commutation relations: 



[fflfcjOfc,] = S k ,k' 



(2) 



where ey& is the Levi-Civita antisymmetric tensor. In 
this model one assumes a cut-off energy A, where A is 
some non-universal quantity that is associated with mi- 
croscopic properties of the bath (A is usually proportional 
to the inverse of the lattice spacing a) . 

The physics described by Hamiltonian (QJ can be sum- 
marize as follows. When S is decoupled from the en- 
vironment (A = 0) one has an isolated spin problem in 
the presence of a "magnetic field" proportional to A. If at 
certain time t = the spin is prepared in an eigenstate of 
Si, the "magnetic field" induces transitions between the 
eigenstates of Si and the expectation value of the op- 
erator Si, namely, (Si(t)), oscillates harmonically with 
frequency A. There is no release mechanism for the en- 
ergy in the spin. By switching on a small coupling to the 
bath of oscillators, the harmonic oscillations of (Si(i)) 
become underdamped due to the dissipation. Second or- 
der perturbation theory indicates that the behavior of the 
system depends on a dimensionless coupling a — X 2 /8tt. 
For a < 1/2 there are two main effects^: the slow modes 
of the bath, that cannot follow the motion of the spin, 
lead to damping and therefore to an exponential decay of 
(Si(i)); the fast modes of the bath, that can follow the 
motion of the spin, lead to a new renormalized oscillation 
frequency Ar < A. For a > 1/2 there is a crossover to 
an overdamped regime where oscillations disappear (ef- 
fectively Ar — * 0) and only exponential decay occurs. 
Finally, at a — 1 there is a true quantum "phase transi- 
tion", where the the impurity spin becomes localized in 
one of the eigenstates of Si . In the Kondo language the 
change from delocalized to localized is equivalent to a 
Kosterlitz-Thouless transition (KT) between the Kondo 
problem with ferromagnetic coupling (that has a triplet 
ground state) and the Kondo problem with antiferromag- 
netic coupling (with a singlet as ground state). 

One of the most illuminating ways to describe the KT 
transition is via a perturbative renormalization group 
(RG) calculation in leading order in A/ A <C 1. The 
RG proceeds in two steps. In the first step one reduces 
the cut-off energy of the bosonic bath from A to A — dA 
by tracing out high energy degrees of freedom. In a 
second step the dimensionless coupling constants a and 
h = A/A are rescaled to the new cut-off leading to the 
RG equations^: 



da 

~di 

dh 
dl 



-h 2 a . 



(l-a)h 



(3a) 
(3b) 



where dl = dA/A. Thus, for a > 1 the system scales 
under the RG to weak coupling (h(£) — > 0), and at low 
energies the tunneling splitting A (£) scales towards zero 
leading to localization. Conversely, for a < 1 the cou- 
plings scales towards strong coupling (h — > oo) indicating 



that RG breaks down. The renormalization scheme fails 
at a certain energy scale (that is, the value of I = £* for 
which h(t*) w 1). This characteristic scale is called the 
Kondo temperature that can be obtained directly from 
Q as: T K w A(A/A) 1 /( 1 - Q ). In the Kondo problem, for 
frequencies and temperatures below Tk there is no remi- 
niscence of the original impurity spin. This is an extreme 
example of decoherence. 

Although the RG equations clearly captures the 
asymptotic behavior of the spin dynamics, in order to 
observe the cross-over from underdamping to overdamp- 
ing, one has to look at the frequency and temperature 
dependence of the spin correlation functions. This is even 
more important in the context of decoherence, since we 
are interested in measuring observables associated with 
the local degrees of freedom, not with the environment. 
In a spin problem, a particular apropos object is the im- 
purity transverse susceptibility that is given by: 



dt 
2T 



f ([Si(t),Si(0)]) 



(4) 



The imaginary part of x(w), ls a measure of the 

amount of energy that is dissipated from the spin into the 
environment. In the absence of coupling to the environ- 
ment (A = in QJ) we have \ (u>) oc 6(u> — A) indicating 
the spin "oscillates" freely with frequency A. When A > 
two different effects occur in the frequency behavior of 
x"(oj)/u>: (1) instead of a Dirac delta function one finds 
a broadened peak and \"{bj)/uj becomes finite at w = 0, 
indicating that the oscillations become damped; (2) the 
maxima moves from A to a renormalized value Ar due 
to "dressing" of the spin by fast environmental modes. In 
the DTLS, the value of x (u)/to\ u — >o and its width Su> are 
set by the Tk- X (w)/w|u^o oc 1/T k and 5ui oc Tk- In 
particular, in the overdamped regime (a > 1/2) the peak 
in x ( w ) at finite frequency vanishes completely leaving 
a smooth function centered around u = Oii. 

In this paper we are going to study a model that can 
be considered a generalization of the DTLS Q: 



H 



* (4°* + b l b k) + AS 3 

^ Vk |AiSi(afe - a\) 



k>0 

i 



k>0 



+ \2S 2 (b k -b{)} 



(5) 



where there are two independent dissipative baths la- 
beled by operators au and bk with couplings Ai and A2. 
Notice that reduces to the DTLS, eq. (QJ, when one 
of the couplings Ai or A2 vanishes. At first sight, the 
only apparent difference between (3J and Q is the ex- 
istence of an additional bosonic bath coupled to a third 
spin component. Thus, naively one would expect an en- 
hancement of decoherence in comparison with the DTLS 
since more heat baths are present. This naive argument 
fails to grasp that both baths are "competing" with each 
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other for the "ordering" of the impurity. While the cou- 
pling Ai "tries" to localize the spin in an eigenstate of 
Si, the coupling A2 also "tries" to localize the spin in 
an eigenstate of S2. However, we see from <|2J that the 
operators Si and S2 do not commute with each other 
and therefore one cannot find a common eigenstate for 
the spin to localize in. This purely quantum mechanical 
effect leads to a less decoherent environment. We will 
show that when Ai,2 = A and A/A <C 1 the spin dynam- 
ics is always in the underdamped regime, regardless of 
the bare value of the coupling constants. In our previous 
publication we called this state of affairs the "quantum 
frustration of decoherence"—. 

The Hamiltonian (0 was originally obtained in the 
study of an spin 1/2 impurity embedded in an environ- 
ment of large spin-S in d = 3 dimensions^. The mapping 
between these two problems is given in appendix EI The 
magnetic environment has two effects in the dynamics of 
the impurity. The molecular fields produced by the envi- 
ronmental spins favor the alignment of the impurity spin 
along the ordering direction giving rise to a "magnetic 
field" proportional to A. The transverse magnetic fluc- 
tuations (spin waves) produce quantum fluctuations that 
tend to misalign the impurity spin leading to couplings 
proportional to Ai t 2 and therefore to dissipation. In an 
ordered antiferromagnetic spin environment the low en- 
ergy, long-wavelength excitations, are two massless Gold- 
stone modes (two transverse magnon excitations) that 
couple to the two different components of the spin as in 
<EJ . The problem of impurities in magnetic media, espe- 
cially in the paramagnetic phase, has received a lot of at- 
tention in the context of quantum phase transitionsi^iii. 
As we are going to show in what follows, the effect of 
quantum frustration occurs at finite energies or frequen- 
cies and therefore before the asymptotic regime is reached 
(very low frequencies) and the impurity spin fully aligns 
with the environmental spins. Thus, as in the case of the 
Kondo problem, quantum frustration is a crossover phe- 
nomenon that cannot be obtained in "asymptopia". We 
should stress, however, that the phenomenon of quantum 
frustration is more general than its origin would imply. 
As in the case of the Kondo effect, it represents a uni- 
versality class of impurity problems where decoherence is 
reduced by pure quantum mechanical effects. 

As mentioned above, impurity problems can be treated 
by powerful theoretical techniques when reduced to one- 
dimensional models with a boundary. It is convenient, 
therefore, to rewrite (3J in a real space representation: 

/+00 
dx (d x Mx)) 2 + AS 3 
■°° 0=1,2 

- v^^0x&(O)£i-v^^0x^(Q)£a, (6) 

where (f>ifl(x) are one dimensional chiral bosonic fields 
(that is, left movers only) associated with the bosonic 
modes a k (b k ) and we have defined cki^ = Af 2 /87r. We 
are ultimately interested in the general problem of de- 
coherence described by © or © and the mechanism of 



quantum frustration associated with this model. 

The paper is organized as follows: we derive the main 
RG equations in Section [H] and show that the dissipa- 
tive model discussed here is always coherent and shows 
scaling at strong coupling; in Section IIIII we study the 
impurity susceptibility using numerical renormalization 
group and analytical RG via the Callan-Symansky equa- 
tions; Section IIVI contains a discussion of the problem 
of frustration of decoherence and also our conclusions. 
There are various appendices where the details of the 
calculations have been included. 



II. RENORMALIZATION GROUP 

Notice that, according to the RG equations (0), the 
KT transition occurs at a finite value of the coupling con- 
stant a and therefore cannot be obtained directly from 
perturbation theory. Instead, one has to use a rotated 
basis of states, obtained from a unitary trasnformation, 
where the problem becomes perturbative. This can be 
accomplished in our case by defining two unitary trans- 
formations: 



JJ _ e i%S 2e i\/8Tra 1 <pi(x=0)S3 
U-2 



(7a) 
(7b) 



that rotate the impurity spin around the S3 direction 
by angles that depend on the field configurations and 
around S2 (Si) by ir/2. Notice that Ui (U2) generates 
a non-perturbative rotation in terms of the coupling ot\ 
(a 2 ). 

Let us consider the problem after rotation by U%, By 
applying U\ to the Hamiltonian ©, we obtain 

ulHUx =H + ^ (AA+ + iy/S^B+) + h.c. (8) 

where Hq is the free bosonic Hamiltonian (the first term 
in the left hand side of ©). We have defined two vertex 
operators, 



_ e =FiV87rai 01 O=0) g± 



Bf = d x <h = 0) e^Vs^M^o) S d 



(9a) 
(9b) 



where S ± = Si ± iS% are the standard raising (lowering) 
operators. 

As in the case of a generalized Coulomb gas 
problemi&»A£, the partition function of the problem, 
Z, can be obtained in the basis that diagonalizes S3 
(Ss\s 3 ) = ±||s 3 )) as 



irrijy/STTa^B'^ 3 (r,)] 



(10) 



where S is the action for the free boson fields, St is 
the time step in the imaginary time direction, and rrij = 
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s 3 ( T j + St) — S3 (tj) is either +1 for a kink or —1 for an 
anti-kink at time Tj of a given spin history in imaginary 
time. The partition function given in iflOj) is the starting 
point of the RG analysis. 

We can define the Fourier transforms of the vertex 
operators, A\{uj) — J dr exp{iuiT}Ai(T) and bosonic 
fields (f>i t 2(k,u)) = J dx J dT(j>i^(x, r) exp{i(kx — lot)}, 
and divide the fields into slow modes, say Ai t< (r), with 
(lu, k) < A and fast modes, say Ai t> (r) with (to, k) > A. 
We then integrate the fast modes within a shell A < 
(k, to) < A + dA, to obtain the renormalization of the 
slow fields due to the fast modes. In this procedure the 
renormalization of the slow modes is given by averages 
over the fast modes. It is straightforward to show that: 



(Af(T)) > = ^We-^ 
(B 1 ± (r)) > = Bf t< (T) 



,-(l+Ql)0 



(11a) 
(lib) 



where, dA/ A = d£ and (P)> indicates the average of the 
operator P over the fast modes. Substituting dill) into 
ltlH)l and rescaling the fields in order to obtain the same 
partition function with slow modes only, we find that 
the couplings have to change with £ according to (see 
Appendix El : 



da 2 

~dJ 
dh 

~dl 



—■10L\OL 2 

(1 - ax)h, 



(12a) 
(12b) 



which define the RG equations for a 2 and h but not for 
ol\. The RG equation for a% is obtained in second or- 
der in h. In the language defined by Anderson- Yuval- 
Hamanni^, it corresponds to the renormalization in a% 
due to a "close pair" of flip and anti-flip that is removed 
from a spin history in a particular RG step. One can show 
that a new operator, which is not present in the original 
problem is generated under this procedure^. This oper- 
ator reads: 



C » 1 - iy/2^h 2 d T (/) 1 (0, t) S 3 (r) diSi 



(13) 



This term can be reexponentiating into the action, 
Eq. ([Tfi)) . and then integrated by parts in r. The final 
result is equivalent to a redefinition of the vertex opera- 
tors, 



At 



o =R\/8^F^r(l-i/i 2 cK)0i(O) g± 



(14a) 



d x <f> 2 (Q)e TtVW ^( 1 -^ h2de )' t ' 1< - 0) S ± , (14b) 



immediately implying the RG equation for a-^, 
da x 2 

~dT ^ ~ ai - 



(15) 



Eqs. ill 211511 where derived by a perturbative treatment 
in powers of a 2 and h and are valid up to second order 
in these coupling with a\ being arbitrary. If instead we 
apply the unitary transformation Eq. l|7bjl a similar set 



of equations can be derived for ot\ and h small with a 2 
being arbitrary. Notice that the only change in the RG 
equations is the interchange between a\ and a 2 in (I12al 
11511 . In fact, given the form of the the Hamiltonian Q it 
is easy to see that the RG equations must be symmetric 
under the interchange of ol\ and a 2 . Thus, it is straight- 
forward to see that by symmetry the RG equations are: 



da>i 

~dl 

da 2 

~dJ 
dh 

~d£ 



—2a\a 2 — aih 2 , 
—2a 2 a\ — a 2 h 2 , 
(1 — ai — a 2 )h. 



(16a) 
(16b) 
(16c) 



The symmetrization process is just a simple way to obtain 
the next order corrections to the RG equations. Strictly 
speaking, the RG equations 11 (ill are valid up to second 
order in h, when either both ai and a 2 are of the same or- 
der and small, or when one of them small and the other is 
arbitrary. However, the terms of the form aia 2 could also 
be directly obtained from a diagrammatic technique^. 
Notice that in the highly anisotropic case, say a 2 = 
(«i =0), we identify ct\ = a (a 2 = a) so that eq. I|16a|) 
(eq. l|16bjl l reduces to eq. 1-iall and eq. II16cH becomes l|3bll . 
As expected, our problem maps into the DTLS and one 
obtains a KT transition at a\ = l (a 2 = 1). The RG 
flow associated with eqs. ltl6)l in the ai versus a 2 plane 
for fixed h is shown in Fig. 




Figure 1: Renormalization group flow given by eqs. in 
the a\ versus a.i plane. 



In the fully symmetric case where a.\ = a 2 = a one 
finds a very different physics. Indeed, from 11 (it . one 
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gets: 



da 

~dl 

dh 

~d£ 



= -2a 2 - ah 2 , 



= (l-2a)h. 



(17a) 
(17b) 



As one can see from Fig. there is no KT transition in 
this case. The couplings a\ and a 2 always flow to zero 
while h scales towards strong coupling. In the DTLS 
language the spin never localizes in an eigenstate of 5i 
or 52 being always in an eigenstate of 53. Hence, in 
the isotropic case, no matter how large the couplings to 
the environment the spin is always coherent. This is the 
phenomenon of quantum frustration of decoherence. 

We can obtain a more quantitative analysis of the RG 
scale in some particular limits. As noticed before the RG 
breaks down at a scale £* = ]xl(Aq/Ta) (where Ao is the 
initial cut-off of the problem) when h(£*) w 1. Ta is 
the crossover energy scale from weak to strong coupling 
(the equivalent of the Kondo temperature) . It is easy to 
see that the value of Ta depends on the bare value of 
a(£ — 0). If a(0) <C h(0) the flow is essentially the same 
as the usual KT flow and one can disregard the flow of 
a(£) in order to find, 



2q(0)/[1-2q(0)] 



(18) 



which is a valid result even when 2a(0) ln(Ao/A) ~ 0(1) 
although its derivation requires a(0) <C 1. If, on the 
other hand, a(0) > h(0) then the a 2 term dominates 
and the flow of a{l) and we must take into account the l- 
dependence of a{l) in solving for the flow of h(l) to strong 
coupling. This leads to: 



Ta ~ A(l + 2a(0) ln(A / A))" 



(19) 



Observe that II18J1 and l|19|) are identical when 
2a(0) ln(Ao/A) <§C 1 but give a very different result when 
2a(0) ln(A /A) ~ 0(1). We immediately notice that the 
a 2 term in the RG destroys the KT transition. Unlike 
the Kondo problem the system retains coherence even at 
large coupling and is never overdamped. This is a quan- 
tum mechanical effect and comes from the fact that the 
spin operators do not commute. While the 5i operator 
in (01 wants to orient the impurity spin in its direction, 
the same happens for the 52 operator. In a classical sys- 
tem (large 5) the spin would orient in a finite angle in 
the XY plane. However for a finite 5 impurity this is not 
possible and the impurity coupling is effectively quantum 
frustrated reducing the effective coupling to the environ- 
ment. Another interesting feature of the RG flow is that 
for h(£) —>■ we find, 



a(t) = 



a(0) 



1 



l + 2a(Q)£* 21n(A /T A )' 



(20) 



when 2a(Q)£* >• 1, a(l) is essentially independent of a(0) 
at energy scale Ta- While Ta gives the crossover energy 



scale between weak and strong coupling, a* provides in- 
formation about the dissipation rate, r _1 , of the impu- 
rity dynamics. Our results indicate that for a(Q)£* suf- 
ficiently large, r _1 is independent of the initial coupling 
to the bosonic baths. 
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Figure 2: Renormalization group flow given by eqs. 11711 in 
the a versus h plane. 

In Fig. we depict the RG flow in the a versus h 
plane. As discussed above, we can see that asymptoti- 
cally (that is, large £) a{£) renormalizes to zero while h{£) 
becomes large. An interesting feature of this RG, as we 
pointed out above, is that for large values of a(0) (large 
coupling to the environment) and intermediate values of £ 
the renormalization of h{£) becomes independent of a(0). 
This indicates that there is a single variable that deter- 
mines the RG flow at intermediate energy scales. The 
fact that only one coupling determines the RG flow indi- 
cates that there must be scaling in the physical properties 
with the renormalized value of h. In the next section we 
will discuss how the RG results reflect on the behavior of 
the transverse susceptibility. 



III. IMPURITY SUSCEPTIBILITY 

In the previous section we discusse the RG calculation 
in the weak coupling limit. The RG indicates that for 
large values of the couplings nothing new should hap- 
pen. Nevertheless, given the perturbative nature of our 
analysis, this conclusion may not be warranted. Our con- 
clusions can be put on firmer ground with the use of nu- 
merical renormalization group (NRG)£. In NRG we do 
not look at the renormalization of the couplings, as we 
did in the previous section, but at the behavior of the 
susceptibility itself. Thus, in the first part of this section 
we study the behavior of the susceptibility as a function 
of the frequency at T = with NRG. In the second part 
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of this section, based on the perturbative RG of the previ- 
ous section and the NRG, we obtain analytic expressions 
for the transverse susceptibility in various regimes. We 
show that these two methods provide full support for the 
RG equations obtained in the previous section. 

A. Numerical Renormalization Group (NRG) 

In order to learn more about the model we have 
performed numerical renormalization group (NRG) 
calculations^ on the Hamiltonian ©. Although NRG 
has recently been extended to bosonic models 20 , we fol- 
low a more traditional approach and transform © into 
a fermionic problem. However, the bosonic baths <j>\ and 
4>2 being Ohmic, we can also represent them as the spin 
density fluctuations of two fermion fields, ipi and ip2, 

A x - t 

k,fj,,i=l,2 

+ 7;.9i<S'i'0l cr iV'i + 7j92 Si(>l a 2 1P2 , (21) 
where vf is the Fermi velocity and Vv = J2 c ikn are 

k 

the local fermion operators. Notice that we have two 
different set of fermions (labeled by i = 1, 2) that couple 
by x and y-component of their "spin" to the corresponding 
components of the impurity spin. 

In order for ll2Dl to be a faithful representation of 
© one has to map the bosonic couplings ai,2 into the 
fermionic couplings gx^. As in the case of the Kondo 
problem^ the bosonic couplings are related to the elec- 
tronic couplings through the electronic phase shifts 8\ 
and 62: 

ai = 2 {^k\ ■ ( 22 ) 

Here the phase shifts can be determined directly from 
the NRG spectrum. The price what one has to pay 
for this simplicity is that the entire parameter space of 
the fermionic model < gt < 00 covers only a smaller 
regime of the original model < a>i < 1 and therefore 
the localization transition is beyond the boundaries of 
the method. The phase shifts are given with a very good 
accuracy by: 

Si = atan(/(A NRG ) 5l ) , (23) 

where Anrg is the parameter of the logarithmic dis- 
cretization used in NRG and /(Anrg) is a numerically 
determinable factor close to unity. For the numerical 
work we used Anrg = 2 and we find /(Anrg = 2) = 1.03 
(see Fig. EJ, 

In Fig0]we show the results for x±( u> )/u (normalized 
to its value at w — 0) as a function of w/T^ (where 
is the crossover energy - see previous section) in the case 
when a\—a.2 — OL (gi — 52) as one varies a. Notice that, 
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Figure 3: The phase shift (and therefore the bosonic coupling) 
extracted from the NRG finite size spectra as a function of 
the fermionic coupling. 

in agreement with the RG calculation, the susceptibility 
retains a peak even for strong coupling indicating that the 
spin remains coherent. Furthermore, as the coupling in- 
creases the susceptibility curves collapse into a universal 
curve showing that at large couplings to the environment 
the susceptibility can be written in a scaling form: 

X>,^)=Xof(^) (24) 

where \o = d^X^i^ = 0, a, h) and f(x) is a universal 
function so that i{x — » 0) = x and {(x — > 00) i=s 1/x. 
These results are in agreement with our earlier conclu- 
sions based on the RG calculation. 



• a t =a z =0. 15 
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Figure 4: x'K^/u as a function of w/Ta- 
To compare results for our model with that of the sin- 
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gle bath DTLS, we have calculated X±( u; ) f° r a i — 0-59 
and «2 = and compared with the case where ol\ = 
a 2 = 0.59. The result is shown in Fig0 Notice that in 
the DTLS case there is no trace of the peak in the suscep- 
tibility indicating that the relaxation of the spin is com- 
pletely overdamped. However, in the isotropic case one 
finds a well defined peak even when the coupling to the 
environment is large, indicating that the spin still keeps 
memory of the tunneling splitting, even when strongly 
interacting with the bath. This is a clear demonstration 
of the effect of frustration of decoherence. 



3.0 



a,=0.59 a„=0 



A/A=0.01 
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0.0000 
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Figure 5: x'± ( LJ )/ <jJ as a function of lu/A. 



B. Analytic Results 

The RG results of section [H] show that the transverse 
couplings of the impurity to the environment always flow 
to a — > indicating that a perturbative approach should 
give a sensible result. When a — the ground state 
of the problem is an eigenstate of S3 and therefore the 
transverse susceptibility has a Dirac delta peak at uj = A, 
that is, zero relaxation rate, r _1 = 0. In order to obtain a 
finite relaxation one makes use of the Bloch equations^ 
for the expectation values of the spin operators, M, = 
(SiY- 

dM A -j _ M l x + M 2 y M 3 z 

= — M x z , 

dt 2 T 2 Ti 

where I/T2 is the transverse and 1/Ti is the longitudinal 
relaxation rates. It is straightforward to write a second 
order differential equation for Mi (t) : 

d 2 Mi 2 dMi /A 2 
-l + I — 

4 



dt 2 



T 2 dt 



1 



— + ^]M 1= 0, 



implying that the transverse correlation function has the 
form: 



x'i M 



2/T 2 



(uj 2 - A 2 /4 - l/T 2 y + Aui 2 /T, 



(25) 



In appendix IC 31 we derive these results using a random 
phase approximation (RPA) and improve on them by re- 
placing the bare values of the parameters by their renor- 
malized RG value: 



X±H = [arctan^r)]- 1 ^^ 
w (oj 2 - (T A ) 2 - 1/t 2 ) 2 + 4w 2 /t 2 



where 



-^(a*) 2 T A . 



(26) 



(27) 



Notice that 1(26)1 reduces to a Dirac delta function at 
u; = A as a(0) — > 0, as expected. We find that this ap- 
proximation is good for w <C Ta and also describes well 
the NRG results for all u < A when A > T A > A a* . 
In the zero frequency limit lj26ll reduces to 



X^ = 0)^(a*) 2 u/(T A ) 2 + O[(a*) 4 



(28) 



and the Kramers-Kronig relation immediately leads to 
real part of the susceptibility: 

x'Aoj^O) = 7r/[8T j4 (l + (a*) 4 )arctan(l/(a*) 2 )] 

w \/{AT A ) + 0[{a*) 2 ]. (29) 

Although the RPA result gives good results in certain 
regimes it fails in the asymptotic cases. In those regimes 
a new approach has to be developed. For that purpose 
we will use the criteria of renormalizability of the theory 
in order to calculate the susceptibility. If we knew the 
exact /3-functions of the theory, 



(30) 



one could in principle integrate the exact RG flow in order 
to obtain the exact result. However, we only have access 
to the perturbative result llfill that indicates that there 
is no other fixed points in the problem. The question is 
whether these results are valid in other regimes. 

Let us consider some limiting cases of the problem at 
hand. Firstly consider the special situation where a (0) = 
«i (0) = a 2 (0) and there is no magnetic field, A = 
(h (0) = 0) . In this case the Hamiltonian of the problem 
can be written, from (JJ, as: 



H efi = I dx (d*Mx)f- V8^d x M0)Si 

J ~ 1 n 



a = l,2 
^d x (j) 2 (0)S2 



(31) 



From the renormalization group equations, Eqs. II17)I . we 
find: 



/3(a) = ^ = -2a 2 + 0[a 3 ]. 



(32) 



At finite temperature T< A the RG flow is cut-off by 
the temperature and we can write dl w —dT/T and use 
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the temperature as the cut-off. We can solve i'.V2l) for a 
as a function of T at once: 



When T 



a(T)« 

one finds: 

a(T) 



1 + 2a ln(Ao/T) 



21n(A /T)' 



(33) 



(34) 



which is independent of a(0) in agreement with ll2(H . 

We first consider the susceptibility at finite T and zero 
frequency, 



x(T,A,ao) 



4T 



g(T/A,a ), 



(35) 



where g(x) is a dimensionless function. Since the the- 
ory is renormalizable, the susceptibility should obey the 
Callan-Symanzik (CS) equation^: 



g(a,T/A) = 0, (36) 



where 7(a) is the anomalous dimension associated with 
the operator o\. Equation (|36l) expresses the fact that 
a change in the cut-off, A, can be exactly compensated 
for by a change in the bare coupling, a, together with a 
rescaling of the susceptibility. The most general solution 
of CD) is: 



g(a,T/A) = cxp 




[2j(a)/P(a)]da 



h[a(T)}, 



(37) 

where h{x) is an arbitrary function of the renormalized 
coupling. We can rewrite I|37J1 in a slightly different form: 



g(T) = exp 



[2 7 (a)//3(a)]da 



[a(T)], (38) 



where we have introduce a new function <f>[a] and used 
that 



exp 



«(T) 



[2-y(a)/(3(a))dc 



is by itself some (in general unknown) function of 
a(T) and we have absorbed this term into the func- 
tion h(a(T)). Hence a non-zero anomalous dimension 
implies some residual explicit dependence of g(T), and 
hence x(T), on the bare coupling ao in addition to its 
implicit dependence on ao through the renormalized cou- 
pling. Notice from 13411 that a(T) becomes small at low 
T and therefore one can expand 4>(a{T)) in a power series 
in a(T). In this case, replacing il38)l into (l35l) we find: 



X{T) = i cxp 



[27(a) //3(a)] dc 



a"(T), 
(39) 



where b n are the coefficients of the expansion of <fi[a]. 

Eq. (l39l) is formally exact. However, one does not know 
the anomalous dimension a priori. One way to go about 
this is to compare the exact result I|39J1 with the pertur- 
bative result obtained in leading order in ao. In appendix 
[Owe show that perturbation theory gives: 

x(T) = ^[l + 2a o ln(T/Ao)+0[a^]. (40) 

Replacing i|34l in (|39l) and comparing with (|4Til) we find 
that 60 = 0, bi — 1 and 



[27(a) / /3{a) da] ss — ln(ao). 



(41) 



The coefficient of 2 in front of the second term in 14(11 is 
crucial. Note that what fixes the definition of a is the 
RG equation, Eq. (HSt - Since f3{a) ~ —2a 2 , we see that: 



7(a) = -a + 0[a 3 



(42) 



is the value of the anomalous dimension in leading order 
in a. Therefore, we have concluded that 



X(T) 



1 



8Ta ln(A /r) 



(43) 



when T — > 0. This result is expected to be true even at 
very low T when aoln(A/T) > 1. Suppose that the bare 
coupling, ao is not small. What can we say from the RG 
in this case? As long as we consider very low T where 
a(T) <C 1 so that <f>{a[T)) ~ a(T), we have: 



X(T) w exp 



[2 7 (a)//3(a) 



1 



8Tm(A /T)' 



(44) 



The first factor is some unknown function of the bare 
coupling but the T dependence is the same as before. 

Now consider the susceptibility at T = but finite 
frequency. Once again, thanks to the renormalizability 
of the theory obeys the same CS equation with 

the same /3-function and the same anomalous dimension, 
7(a). This anomalous dimension is a property of the spin 
operator o\ and must be the same for either finite T and 
oj = or finite 10 and T = 0. Therefore, following the 
earlier discussion it must have the form: 



x"(w, A, a Q ) = exp 



[2 7 (a)//3(a) 



1 



F(aM). (45) 



The function F(a(ui)) is not necessarily the same as 
4>{a(T)) and, in general, is unknown. However, the first 
factor, giving the explicit dependence on ao in a pertur- 
bative expansion, should be exactly the same as in the 
previous calculation of a(T). Thus, if ao <C 1, we must 
have: 



x"{u,A,a ) = —F{a(u)). 
ujao 



(46) 
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Again, we perform ordinary perturbation theory for 
X (w) in powers of an and improve the perturbative re- 
sult with the RG by matching it to lUfil) by expanding in 
powers of a(uj). Since we already know (3(a) and 7(a), 
the result must have a rather restricted form to be a so- 
lution of the CS equation. The susceptibility at finite 
frequency is given by Eq. itC5j) : 

X "(lo, A, ao) w - [ao - 4a 2 In(Ao/w) + C(ag)] . (47) 

This result is consistent with the RG form of Eq. ijlfif if 
we assume that: 



F(a(u)) w a 2 (w), 



(48) 



On 



Aal m(A /w) 



Having found the function F(a(u>)) at small a(w) we can 
now invoke the RG. In particular, for small bare coupling 
and small w we have: 



1 



4anwln (Aq/lv) 



(49) 



Even if the bare coupling is not small, but we go to small 
enough u so that a(uj) -C 1, the RG implies that: 



x"(uj) » exp 



27(a) / /3(a)da 



1 



4wln 2 (A / 



(50) 



where the first term in an unknown function of the bare 
coupling constant. Thus, eqs. and lH9|l give the tem- 
perature and frequency behavior of the susceptibility for, 
a* <C 1, Ta <C Ana*, and in the frequency and tempera- 
ture range Ta <C oj,T <C Ao. When these conditions are 
satisfied the ratio, 



X '>,T = 0) _ 2T ln(A /T) 



X(T) 



u ln 2 (A /w) ' 



(51) 



is universal. 



IV. DISCUSSION AND CONCLUSIONS 

Decoherence can be defined as the unavoidable evolu- 
tion of the total state of the system and the environment 
towards an entangled state. This is a dynamical defi- 
nition of decoherence, and clearly shows the conceptual 
difference between dissipation (that involves the transfer 
of energy from the subsystem to the environment) and 
decoherence. An important concept in the study of de- 
coherence is the notion of a preferred basis: every time 
that a system interacts with an environment, a set of 
states is naturally selected by the form of the interaction. 
The text book example is given by the exactly solvable 
model22i2i 



Hdeco = ka\ak + % 



fe>0 



■^=Si5^V*(a*-4). 
^ 2L fc>0 



Although this model does not have any dissipative mech- 
anism, the two level system experience strong decoher- 
ence. Suppose that the system is prepared at time t = 
as a direct product of the bath and the two level system: 

p{t = 0) = Pbath(t = 0)® P§ (t = 0), 

where pbath and pg are respectively the density matrices 
of bath and the two level system. A natural basis choice 
for the two level system is Si. If we further suppose that 
the two level system is prepared in a state of S2 , then at 
t = the reduced density matrix, 

Pn(t = 0) = tr bath [p b ath(t = 0) <g> pg(t = 0)] , 

has off-diagonal matrix elements indicating that the sys- 
tem is coherent. As the system evolves in time the off- 
diagonal elements decay very fast due to the entangle- 
ment of S and the bath degrees of freedom. As t — > 00 
only the diagonal elements of the reduced density ma- 
trix of the system remain, and we say that the system 
"decoheres" to the preferred basis of Si . 

With this example in mind is now simple to under- 
stand the effects that we described in this manuscript. 
Consider the Hamiltonian (0 . In this case it is no longer 
possible to define a preferred basis for the two level sys- 
tem. The entanglement of S with each one of the baths is 
suppressed by the other, and as a result the decoherence 
phenomena is frustrated. This physical picture shows the 
true meaning of our results, the "quantum frustration" is 
the lack of a preferred basis for the system of interest. 

The quantum frustration of decoherence can be also 
understood as a result of a version of Coleman's or 
Mermin- Wagner's theorem^. When ai = a2 there is a 
U(l) symmetry in impurity problem. Hence, one has an 
effective (1 + 1) dimensional field theory with U(l) sym- 
metry so this symmetry cannot be spontaneously broken 
even at T = 0. In fact, because one has a single bound- 
ary degree of freedom, one can also think of the problem 
as an almost (0 + 1) dimensional field theory. So it is a 
rather remarkable fact that even the Z2 symmetry which 
remains when ai = can be spontaneously broken, as in 
the case of the DTLS. The U(l) symmetry would have 
to be spontaneously broken in a phase in which the spin 
is localized in an eigenstate of either Si or S2. Quantum 
mechanics prevents that from happening. 

One can ask how generic this result really is. For quan- 
tum frustration to occur the coupling constants with all 
the baths must be identical. This can be achieved when 
the role of the two baths is played by two Goldstone 
modes, resulting from the spontaneous breaking of a con- 
tinuous symmetry, such that the residual unbroken sym- 
metry rotates the two Goldstone modes into each other. 
When the couplings are not exactly equal quantum frus- 
tration occurs up to a certain energy scale below which 
one of the heat baths takes over and one obtains the stan- 
dard decoherence problem in dissipative ohmic systems. 
In terms of FigQ] it means that the asymptotic flow is 
the one for either ai = or a 2 = 0. In summary, quan- 
tum frustration of decoherence is a general phenomena. 
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It has clear implications to quantum/classical transition 
and measure theory. Moreover it is potentially important 
to the development of technologies where decoherence is 
a fundamental issue as in the case of quantum communi- 
cation and quantum computation. 

In summary, we have studied a model of quantum 
frustration of decoherence in open systems. Contrary 
to standard dissipative models with ohmic dissipation, 
the non-commutative nature of spin operators lead to a 
frustration of decoherence. We have shown that while 
in a DTLS the spin dynamics becomes overdamped at 
large couplings with a heat bath, in a system with quan- 
tum frustration it is always underdamped and the system 
keeps the memory of its quantum nature. Using pertur- 
bative RG calculations we have shown that at large cou- 
plings with the bath the transverse spin susceptibility 
shows scaling with a characteristic energy scale Ta, the 
analogous of the Kondo temperature in the DTLS, that 
separates the region of strong to weak coupling. We have 
supported our claims with NRG calculations and have 
calculated the frequency and temperature dependence of 
the transverse susceptibility using the renormalizability 
of the theory. Our results may be applicable to a large 
class of problems where decoherence plays a fundamental 
role. 
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Appendix A: IMPURITY SPIN IN A 
MAGNETICALLY ORDERED ENVIRONMENT 



In this appendix we will show how quantum frustration 
can arise in the context of a magnetic impurity in a mag- 
netic environment. Let us consider a magnetic environ- 
ment describe by the quantum Heisenberg Hamiltonian 
in the presence of an impurity: 



H 



+ A S ■ s , 



(Al) 



where J is the magnetic exchange between nearest neigh- 
bor spins Si located on a lattice site R4 in d dimensions, 
A is the coupling between the environmental spins and an 
impurity spin S located at the origin of the coordinate 
system. In what follows we will consider the antiferro- 
magnetic case of J > although the ferromagnetic case 
(J < 0) can be studied in an analogous way. 

The partition function of the problem in spin coherent 



state path integral can be written as^: 

Z = J DN 8{N 2 -I) J DrlS(n 2 - 1) e -iS,(N)-S(N,H) 

where N represents the impurity spin and n(r) the envi- 
ronmental spin field, 5 B is the Berry's phase, 



S 



drdr { — 
2.9 



+ \5(r)n(T,r)-N(T)X, 



(A2) 



is the action of the problem where g — c 2 / p s is the cou- 
pling constant (c = 2^/dJas is the spin wave- velocity and 
p s = Js 2 a 2 ~ d is the spin stiffness, a is the lattice spacing 
and s is the value of the environmental spin) . 

Assume that the 0(3) symmetry of the model is spon- 
taneously broken so that the field ft orders. In this case 
we can write: 



n(r,r) ps (ip x (r, r) , ip 2 (r, r) , 1) . 



(A3) 



where cpi t 2 are small fluctuating fields corresponding to 
the two Goldstone modes of the antiferromagnet. A pos- 
sibility that is not considered in this work is associated 
with the formation of a spin texture around the impu- 
rity spin. In a classical spin system a spin texture can 
be formed in the bulk spins due to the presence of strong 
and/or anisotropic interactions. The spin texture can fol- 
low the impurity as it tunnels invalidating the methods 
used here (an instanton calculation is required to take 
into account the collective nature of the texture). The 
results in this appendix are only valid if no spin texture 
is formed around the magnetic impurity. 

In the ordered phase the Berry's phase term is unim- 
portant and can be dropped. Using Eq. l|A3ll the action 
<|A2f) reads, 



S 



drdr < — 
I 2.9 



(d T ^ Q (T,r)) 2 + c 2 (V^ Q (T,ry 



- E 

a=l,2 

+ XS (r) Va ( T , r) ■ N a (r) + XN 3 (r)} . 



(A4) 



We see that the action for the fields (£1,2 is quadratic and 
therefore these fields can be traced out of the problem 
exactly. In this case the effective action for the impurity 
spin becomes in Fourier space: 



9_^_ 
2 



E 

a=l,2 



diodk 



N a (u)N a (-u) 



c 2 k 2 



+ X / drN 3 (r) . 



As we should expect from the spherical symmetry of the 
problem, the angular dependence in k can be integrated 
and we finally obtain 



r(d-i) 



+00 



dq 



d-i 



+ A / drN 3 (r) , 



N a (u)N a (-u) 
(A5) 
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where q = ck, and 



gT(d-l)\ 2 



2 d+2 7r (d-2)/2 r ( d / 2 ) C d ' 



where T(a;) is a gamma function. Integrating IIA5J1 over g 
and Fourier transforming back the frequencies to imagi- 
nary time we find: 



S e ff ~ a 2 J]] / dr / <ir 

a=l,2^ ^ 

+ A y drN 3 (r) , 



rt\d-l 



(A6) 



which shows that the impurity interact in imaginary time 
through a long-range interaction that decays like I/t^ 1 . 

The action l|A6t can be simplified by introducing a 
Hubbard-Stratanovich field that splits the interaction 
term. This can be done with the introduction of one- 
dimensional bosonic fields defined as: 



<j> a (x, t) 



mi — 7= — 4>* a {Q^v 



2L . „ 

g — iqx-\-iuj n T 



4> a (q,w n ) , 



L — > oo is the size of the one-dimensional line. Using 
these new fields the action IjAfiJl can be written as: 

Sen = ^2 T^2^2[iuj n + q}(f>* a (q,uj n )(t) a (q,-uj n ) 

a=l,2 q>0 uj„ 
d-l 

+ a \q\ 2 [<t>*a (q,u n )N a (u n ) + <f> a (g,u„)JV a (-u n )] 
+ A JdrN 3 {r). (A7) 

It is easy to see that the trace over the bosonic fields 
reproduces <l AB|) . It is straightforward to see that in d = 3 
the above action reduces to JfjJ. 



Appendix B: RG EQUATIONS 

In this appendix we will derive the RG equations ll"12t . 
From equation l|9ajl we have: 



(Af} > = ( e = F i V / 8Tai(0i,<(z=O)+0i, > (a;=O)) 
_ j{± ^is/Snai <f>i x> (x=0) ^ ^ 

= A^ < e' aide ^ A^il - ai d£) . (Bl) 
Substituting (IB 1 1) in ill( )fl one obtains a term of the form: 

STe de Ae' aiM A"% = 5tA(£ + d£)A^ 3 < , (B2) 

where we have used that, by rescaling u — > ui/b (with 
b — e de w 1 + di) one has r — ► br. Hence, 



and defining the dimensionless coupling, h(£) = A(£)/A, 
one obtains l|12hj) . 

Analogously, from Il9hjl we have: 



= Bf e -(1+Q!l)< " 

« B^Il-Cl + ax)^], 



(B4) 



where we have used that d x —> b 1 d x since — » fc/fe. 
Replacing (IB4tl into the second term in the r.h.s. of IKH : 



JrV^e-" 1 ^^ ( Ti ) - (Sr^aa^ + dQB^ , (B5) 
and hence we write: 

o; 2 (| + d£) = a 2 (£)e' 2aiM m Q 2 (€)(l - 2 ai d£) , (B6) 
leading to equation Ill2all . 

Appendix C: PERTURBATION THEORY 

In this appendix we show how to derive the perturba- 
tive expansion for the transverse susceptibility: 

S (r) = (T T Si (r) 5! (0)) . 



1. Static susceptibility u = and h = 

Firstly, let us consider the case of arbitrary a 2 but 
small a.\ (the case of arbitrary a\ and a2 -C 1 is com- 
pletely analogous) . This regime can be obtained by using 
eq. 0: 

i?2 = U 2 1 H ei fU2 

= H a -V^a 1 d x <f> l e- l " /W ^ 2 ' t > 2{0) S + +h.c. 
In this rotated basis, S (t) has a simple form: 

S{t) = \(T r (AUr) + A^(r)){A+(0) + A^(0))), 
where 



The leading order terms in an expansion in powers of a\ 
at T = can be immediately obtained from the bosonic 
propagator: 

S{t) « ilDrl-^+Otaxaa], (CI) 

where I? is a short time cut-off. We can use the stan- 
dard conformal transformation to promote this result to 
a finite temperature expression, 



A(£ + di) = A{£)[1 + (1 - ai )d£] 



(B3) 



S(t) rs ~ ( — sinkTr 
W 4 \ Tvr 1 



-2a 2 



C[ai«2] 
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Expanding the above expression for 02 < 1 gives: 



S(r) 



1 



I - 2n 2 ln ( ^sin|7rrT| ) +0[a l a 2 ,a 2 2 ] 



To this order, the susceptibility can be calculated imme- 
diately: 



X (T) = drS{r) 
Jo 



1 

4T 



For completeness, let us re-obtain this result by a direct 
perturbative calculation in second order in a 2 : 



(C2) 



,1/T ,T 

S (r) w 8na 2 J dr\ J dr 2 [(T T a 2 (n) cti (t) a 2 (r 2 ) 01 

l/T 



7ra 2 



(0)) - (T T <7 1 (r) 01 (0)) (T r a 2 (n) a 2 (r 2 ))] (T T d</> 2 (n) (r 2 )) 



dn / dr 2 (50 2 (n) d(j) 2 (r 2 )) 
Jo 



Using the finite temperature propagator, 

1 1 



(n)^ 2 (r 2 )) = 



we obtain 
5(r) « - 



a 2 
T 



l/T 



4vr [_^8in[7rr|r 2 -Ti| 



1 



dn / dr 2 — 

r+i Jo [^BintTrTlra-nl] 



o 2 

T 



In ( ^- ) - In (sin ,t7Y| i 



G[ai(X2, al] , 



in agreement with jC2L 



lation function in fourth order in the coupling constants. 
From Eq. ijClt we already know part of the result, 



S (r) = ~ - ^a 2 In \Dr\ + \p\ In 2 |Dr| + 0[a x a 2> t*K|C3) 



2. Dynamic susceptibility uj 7^ at T = and h — 

The finite frequency calculation is a little more tedious 
than the previous one. We would like to obtain the corre- 



The remaining contribution to the correlation function 
is a term proportional to a\a 2 . A convenient way to 
derive this contribution is to use Eq. lT7a| and compute 
the result to all orders in a\ but for a2 < 1. In second 
order in a 2 we need to calculate: 



J 



SS(t) = i + 27ra 2 J dT 1 dT 2 (T T S 3 (T)B+(r 1 )B^(r 2 )S 3 (Q))-7Ta 2 J dndr 2 (T T B+ ( n ) 5f (r 2 )> . 



From this point, it is straightforward to obtain the corre- 
lation function and the susceptibility at finite frequency: 



SS(r) 



$X fan ^ 0) 



a 2 



1 

4 ~ 4£) 2 Q1 [! + 2a{\ ai 
1 1 



M \2ai 
I D I 



| 2 Q1 



[1 + 2ai] 2ai J D 2 «i \cu n \ 
{r(l - 2ai)sin [ttoi]} 



1-2c*i 



Expanding for ai <C 1 and ai ln(Z?/|w„|) <C 1 we find: 
naia 2 



5x fan ^ 0) 



(C - 1) + In 



(C4) 



where C « 0.57772 is the Euler-Gamma constant. The 
susceptibility in fourth order in the coupling constants 
is the sum of Eq. IIC4JI and the Fourier transform of 
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Eq. EH, 



X( iuJ n) = 7 S(uJ n ) + 



7T«2 1 



1 — 2 («2 + Oil) In 



2Ca 2 + 2 (C - 1) ai + 0[<4,a{] 



(C5) 



3. Asymptotic regime of h — > oo 

We represent the spin variables in l(3"T)) in terms of two 
spinless fermions: 

51 = (a% + tfa) /2 

5 2 = -i (a f 6- 6 f a) /2 

5 3 = (a f a - 6 f 6) /2 

and add to the action an imaginary chemical potential, 
iujo = iirT/2. Working with this formalism we can use 
Wick's theorem and the standard diagrammatic tech- 
nique. The action is rewritten as 

S = So(<fK,<h) + J dTa*(T){d T -z]a(T) + 

+ b* (r) [d T + z*] b (r) + V27 (0, r) 

- i^/a^d x (j>2 (0,r)]a* (r)&(r) +h.c. , 

where z = iwo + §. We define the following propagators, 



C/^ 0) (iu/„) = 



(iu„ ^ 0) 



lLU n — Z 
-1 



M arctan( ^ 

27T V CJ„ 



for the fields a(r), 6(r), and the boundary field 0i,2(O, r). 

From the propagators we immediately derive the ze- 
roth order part of the susceptibility: 

XQ (iw„) = ^i 0) (^n + *Pn) ^b° } On) 

+Gl 0) (iUn + iPn)Gi 0) (ipn) 

1 



iuj n + h iuj n — h 
h 



2u 2 +h 2 ' 



A simple perturbative calculation will fail to capture 
the physics and the correct behavior of the susceptibility. 
Following the standard prescription we will sum the in- 
finite series of bubble diagrams in the RPA. Let us first 
consider the second order bubble diagrams. From the 



definition of the propagators and assuming |w n | <C D, 
we obtain: 



8xi fan ^ 0) 



5x2 (iu n ^ 0) 



TTttl h . . 

- — sWr 



= -2nai [xo ( 



7ra 2 w„ 



OJr, 



2 ( u * + h*y 

= 2ttq!2 [xo (iu n )] 2 (^-J \uj n \. 

The bubble diagrams in fourth and sixth order can be 
calculated in a straightforward way: 



(4) /■ x 
Xrpa \ luJ n) . 2 2 

~ 47T UJ„ 



[Xo {iu n )f 
(6) n \ 

X R PA i^n) 
[XO {iUn)f 



al~a 2 (^j^) (2ai+a 2 ) 



-8^ 3 |w„| 3 



2 {u n y 



2«1«2 + Ctn) + Oi 



For h ^> a 2 w„, we can simplify these results and sum the 
geometric series, 



Xrpa (iu n ) 



Xo [iu n ) 



l + 2 7 ra 1 |w„|x (0) (»w„)' 

W2) 
h 2 + u 2 + nhai \w n \ 



The zero temperature susceptibility in the RPA approx- 
imation (for low frequencies and high magnetic fields) is 
obtained by the analytical continuation, 



Xrpa M 



a\h 



2 {lo 2 - h 2 ) 2 + n 2 ha 2 uj 2 



If we define the decoherence time 



T,- 1 = -hal 
2 2 1 



(C6) 



we can identify the functional form obtained in Eq. 1(251) , 

XrW H „ h/T 2 

w ~ (w 2 -/i 2 ) 2 + 4w 2 /T 2 2 ' 
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that leads to 12611 if one replaces: 

h - T A 

U\ —> a* 

T,- 1 - T~ x = ^T A {a*f. (C7) 
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